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Abstract 

We find the general behaviour of homogeneous and isotropic cosmological models in 
some fourth-order theories of gravity. Explicit, exact, general solutions are given for both 
empty universes and those filled with a perfect fluid. For the vacuum case, solutions 
are found with closed, open and flat geometries, whilst the perfect fluid solutions are all 
spatially flat. Both early and late-time limits are studied, and attractor behaviour towards 
simple power-law expansion is identified. Multiple solutions to the same theories, with 
the same matter content and topology are found. It is shown that these solutions exhibit 
great variety in their evolution. 
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1 Introduction 



Generalisations of the usual Einstein-Hilbert action of general relativity (GR) have been ex- 
tensively studied in the literature. One frequently considered modification is to replace the 
Ricci scalar, R, by some analytic function, f(R) (see e.g [1, 2, 3, 4]). Theories of this type are 
referred to as fourth-order, as the field equations generated from them are generically fourth- 
order in derivatives of the metric. Motivations for such studies are found from various different 
sources. One often cited reason is that early attempts to create a perturbatively renormalizable 
quantum field theory of gravity found success by adding extra terms, quadratic in the Ricci 
curvature, to the action [5]. More recently, the effective actions of some string theories have 
been shown to include higher-order curvature terms [6, 7]. Studies of f(R) theories have also 
been performed in cosmological settings, often in attempts to better understand the late-time 
accelerating expansion of the universe [8, 9], cosmological inflation [10, 11, 12] or the nature of 
an initial singularity [13, 14, 15]. 

Difficulty can arise, in studies of fourth-order theories, as the field equations involved are 
considerably more complex than their counterparts in GR. This extra complexity brings new 
and interesting behaviour, such as violations of Birkhoff 's theorem [16] or violations of the no 
hair theorems of de Sitter space, [12]. However, this extra complexity also makes it even more 
difficult to find exact solutions of the field equations. Whilst some exact solutions have been 
found in the spherically symmetric [16, 17] and cosmological [17, 18, 19, 20, 21, 25] situations, 
these are all particular solutions. To date, there have been no general solutions published in 
the literature, for any set of non-maximal symmetries. Exact solutions, and particularly exact 
general solutions, are of great importance for understanding a theory. In this paper we will 
present homogeneous and isotropic exact general solutions, obtained through direct integration 
of the field equations, to some fourth-order theories of gravity. We expect these solutions to 
be of use for helping to understand fourth-order theories of gravity, and the evolution of the 
universes they govern. 

Investigations of fourth-order theories usually follow one of two approaches: Either they 
take the full-theory and look for approximate solutions, or they approximate the theory and 
look for exact solutions. We will take the later approach and consider theories that have a 
Lagrangian proportional to R n . Such theories are scale invariant, and reduce to GR in the 
limit n — > 1. These theories may be considered as the limit of a more general Lagrangian 
that has a power of R dominating in some particular regime, or as a simple deviation from the 
standard theory in their own right. We find that for homogeneous and isotropic cosmologies 
the field equations of these theories can often be integrated directly. For spatially flat vacuum 
cosmologies the general solution can be found for any n, for spatially curved vacuum cosmologies 
we find the general solutions when n = 3/2, and for spatially flat perfect fluid cosmologies (with 
equation of state p = (7 — l)p) we find the general solutions when n = 37/(37 — 1) and when 
n = (10 — 37)/(2(7 — 37)). These dependences are shown graphically in figure 1, below. 

Investigations of homogeneous and isotropic cosmologies in R n theories have been performed 
before in [17, 21]. In these papers the authors find power-law exact, particular solutions and 
perform a dynamical systems analysis of the phase space of more general solutions. These 
studies show that the general solutions are often attracted to simple power-law solutions, at 
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Figure 1: A plot of n = 37/(37 — 1) (solid line) and n = (10 — 37)/(2(7 — 37)) (dashed line). 
For these values of n the field equations for homogeneous and isotropic perfect fluid cosmologies 
can be integrated directly. 

both late and early times. We will confirm this behaviour here, and add to the previous work 
by finding explicit expressions for the general evolution of these universes. 

We will now comment on the extent to which the solutions found in this paper can be 
considered candidates for the description of our universe. It has been shown in [16] and [17] that 
gravitational theories derived from a Lagrangian of the form R n are constrained by observations 
in our local universe to have n very close to 1. It has also been suggested (see e.g. [22], [23]) that 
it may not be possible for general f(R) theories of gravity to transition between eras dominated 
by different fluids in the same way that occurs in the usual general relativistic description. 
These results suggest that our neighbourhood of space-time should be well described by a 
theory that is very close to the Einstein-Hilbert one. However, with these constraints in mind, 
there are still good reasons to be interested in deviations of the R n form. Firstly, even if n is 
very tightly constrained by local observations, any small deviation from n = 1 may produce 
significant deviations at early times in the universe's history. This will be shown explicitly 
below. Secondly, theories of the type R + aR n are often considered in the literature. If n > 1 
then we may expect the effective gravitational theory in the early universe to be well described 
by a Lagrangian of the form R n . In such cases a more appropriate description of the evolution 
of the early universe would be given by an R n theory, even though observations from our local 
universe show strong agreement with an effective theory of the Einstein-Hilbert type. 

This paper will proceed as follows. In section 2 we will give the field equations for the 
theory, and show how they can be simplified by recasting them in terms of new variables and 
transforming time coordinates. Simple power-law exact, particular solutions to these equations 
are given. In section 3 we show how the field equations for vacuum cosmologies can be decoupled 
by a further transformation of variables. These equations are then integrated directly, and the 
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solutions given explicitly in terms of the metric. Section 4 follows a similar prescription, this 
time solving the field equations in the presence of a perfect fluid. In section 5 we perform a 
brief analysis of the vacuum solutions, and in section 6 we analyse the perfect fluid solutions. 
Section 7 provides a closing discussion. For readers interested in the solutions, but not the 
derivations, the solutions in sections 3 and 4 are boxed to make them easily identifiable. 



2 Field equations 

We consider here a gravitational theory derived from the Lagrangian density 

C G = -V=gR 1+s , (i) 

x 

where 5 is a real number and x is a constant. The limit 5^0 gives the Einstein-Hilbert 
Lagrangian of GR, and we are interested in isotropic and homogeneous cosmological solutions 
with 5^0. Theories of this kind are dynamically equivalent to scalar-tensor theories (see e.g. 
[24]). 

We denote the matter action as S m and ignore the boundary term. Extremizing 



S = J C G d x + S m , 
with respect to the metric g ab , then gives [19] 

5(1 - 5 2 )R S ^± - 6(1 + 6)R S ^ + (1 + 5)R 5 R ab - l -g ab RR & 

- g ab 6(l - 6 2 )R S ?^ + 6(1 + S)g ab R 5 ^ = |r a& , (2) 

where T ab is the energy-momentum tensor of the matter fields, and is defined in terms of S m 
and g ab in the usual way. We take the quantity R 5 to be the positive real root of \R\. 

We are concerned with idealised homogeneous and isotropic space-times described by the 
Friedmann-Robertson- Walker metric with spatial curvature parameter k: 

rfs 2 = _ dt 2 + fl2(t) f dr 2 + rHQ2 + r2 s . n2 2 \ 
\(1 — nr z ) J 

Substituting this metric ansatz into the field equations (2), and assuming the universe to be 
filled with a perfect fluid of pressure p and density p, gives the generalised version of the 
Friedmann equations 

(1 - 5)R 1+S + 35(1 + 5)R 5 (J + 3~ j - 35(1 - 5 2 )R^ = \(p - 3p) (4) 
-3^(1 + 5)R S + ^ + 36(1 + S)\llf = \ P (5) 
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where, as usual, 

• 2 

„ (X d nj , x 

R = Q- + Q- + Q—. 6 

a a 2 a 2 

It can be seen that in the limit 5^0 these equations reduce to the standard Friedmann 
equations of GR. A study of the vacuum solutions to these equations has been made by Schmidt, 
see the review [20], and a study of the perfect-fluid evolution has been made by Carloni et al 
[21]. In [17] Clifton and Barrow investigated solutions when k — 0, and showed that attractor 
power-law solutions exist. Various conclusions are also immediate from the general analysis 
of f(R) Lagrangians made in [3] by specialising them to the case / = R 1+s . These previous 
studies have focussed on the existence of particular exact solutions, or qualitative investigations 
of the phase space of general solutions. In what follows we shall be interested in determining 
the general evolution of a in terms of explicit, exact solutions. This should lead to fresh insights 
into the evolution of universes described by these theories. 

Assuming a perfect-fluid equation of state of the form p — (7 — l)p, the usual conservation 
equation gives p = poa~ 31 . Substituting this into equations (4) and (5), with k — 0, gives the 
particular power-law exact Friedmann solution, for 7 7^ 0, 

2(1+6) 

a(t)=t~^~. (7) 

Alternatively, if 7 = 0, there exists the de Sitter solution a(t) = e mt . These power-law solutions 
reduce to the usual general relativistic solutions in the limit 5 — > 0. Furthermore, the power-law 
solutions 

5(1+26) 

a(t) = t a-*) (8) 

and 

a(t) = t* (9) 

also exist. These solutions are independent of 7, and therefore of the matter content of the 
universe. The 6 = limit of (8) clearly corresponds to Minkowski space. These three particular 
solutions were both shown in [17] to be attractors of the spatially flat general solution, under 
certain conditions. In what follows we will be able to show this attractor explicitly in terms of 
exact general solutions. 

The equations (4), (5) and (6) can be cast in a simpler form by defining the new variables 
= V3 In R s and a = e^^a. 

Transforming to the time coordinate 

d V = e^dt 

and setting x = 1 we then have two second-order evolution equations, for the two variables a 
and 0, 

6 ^ = _02 y x (2-37)po e v5 (7 -4/3W2 (1Q) 

+ 3^0 = \V e~ X(p - v / 3(7 - V3)po e V3( T -4/3)^/2 (n) 
a 2a J T 
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and the constraint equation 



where 



6% = ^ + V e-» - 4 + J^e^-^l 2 (12) 
a 2 2 or cr 7 

Vq = . - sign(i?) and A — 



(1 + 6) & v ' ^35 ' 

Over-dots here denote differentiation with respect to rj. In rescaling time in this way we must 
remember that the more physically significant cosmological time coordinate t may diverge at 
finite values of the new time coordinate. In such cases we will consider t as the more physically 
meaningful time, and allow the new coordinate to take values only over an appropriate range. 
This set of equations is similar to those obtained in GR, when considering homogeneous and 
isotropic cosmologies with a scalar field in an exponential potential, and a perfect fluid. This 
similarity is due to the conformal equivalence between f(R) theories of gravity and GR (see e.g. 
[25, 26, 27]). One should notice, however, that the terms corresponding to the perfect fluid have 
a more complicated form than may have otherwise been obtained in GR (except for the special 
case of black-body radiation, 7 = 4/3, which is conformally invariant). The corresponding set 
of equations in GR has been solved in vacuum by Russo [28] , and in the presence of a perfect 
fluid by Dehnen, Gavrilov and Melnikov 



3 Vacuum solutions 
3.1 Spatially flat solutions 

In the absence of any matter fields the equations (10) and (11) are identical to homogeneous 
and isotropic general relativistic cosmologies, with a scalar field in an exponential potential. 
The spatially flat general solutions to these equations have already been found by Russo [28]. 
Here we will briefly reiterate the method of solving these equations, and use the results to find 
the general solution to spatially flat vacuum cosmologies in R n gravity. 

Firstly, making the transformation of variables {a, 0} — > {u, v} by the definitions 
a 3 = e v+u and = ^(v - u) 



and defining the new time coordinate 



di] = e^dr 



allows the evolution equations (10) and (11) to be written as 

A \ . 2 W (. A 



ii+|1 -^r-fl 1 + ^l=° 

A V-2 3Vb ( ^ A 
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where over-dots now denote differentiation with respect to r. In terms of these new variables, 
the constraint equation (12) reads 

(15) 



uv 
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We will first solve for two special cases, and then give the solutions for more general cases. 
3.1.1 S = -1/2 

The special case A = a/3 (corresponding to 5 = —1/2) gives the solutions to equations (13) and 
(14), under the constraint (15), as 



u = a + - ln(r + c 2 ) 



and 



v = c 3 + — (r + c 2 ) 



where the q are constants. Transforming these results back, and absorbing constants into 
redefinitions of the coordinates, gives the metric 



ds 2 



-2/3 



dr 2 + r 2 / 3 rfx 2 



(16) 



where dx 2 is the line-element of flat three dimensional Euclidean space. The ± sign here is due 
to the dependence of V on the sign of the Ricci scalar. 



3.1.2 6 = 1/4 

The case A = — \/3 (corresponding to 5 = 1/4) is also special and gives the solutions 



3Vo , so 
u = c A + —\ r + c 5 ) 



and 



v = c 6 + - ln(r + c 5 ) 



i.e. the same as the solutions in the previous case, under the exchange of u and v. Now, 
transforming variables back to the originals and absorbing constants into coordinate freedoms, 
we have 



ds 2 = _ e ±2r» r 2/3 dr 2 + e ±r» dx 2_ 



(17) 



Here the two ± signs must be chosen together. 
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3.1.3 6^ -1/2 or 1/4 



Under the conditions V^A 2 — 3) > and 5 ^ — 1/2 or 1/4, the solutions to equations (13), (14) 
and (15) are 



u — C7 H t= in 



V = Cg + 



V3 



In 



Transforming back to the original variables, and again absorbing constants, gives 



, 2 COS 4 *- 1 T , 2 , o 

(is = 2 d-r + cos 4,5-1 rdx . 

Si n 27+T r 



In the limit 5 — > this solution can be seen to approach Minkowski space. 
Alternatively, when V^A 2 — 3) < the solutions are 



V3 

u = c io H = In 



V = C12 + 



v/3 + A 
V3 



cue 



■('-*> 



In 



'(-*)• 



+ cue 



v^- A 

which in terms of the original quantities corresponds to the metric 



ds2 = Je^c n e-^ 
(e T + c 11 e~ T ) 2S + 1 



dr 2 + (e r - c u e- T ) 



45 

45 " 1 rfx 2 



(18) 



(19) 



Here we have absorbed the constants c w and c 12 into coordinate redefinitions whilst retaining 
c\\. It should be noticed that whilst the magnitude of en could have been absorbed, it has 
been left as the special case en = corresponds to the attractor solution. Again, this solution 
approaches Minkowski space as 5 — > 0. 

The conditions imposed upon Vq and A here do not imply any restriction on 5. This is not 
immediately obvious as the sign of V depends on the sign of R, and hence on the solution. 
However, calculating the Ricci scalar for solutions (18) and (19), and substituting this into the 
definition of V , shows that both of these solutions exist for any given 5 —1/2 or 1/4). 
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3.2 Spatially curved solutions 
3.2.1 8 = 1/2 

We will now present the general solution for a spatially curved vacuum universe, when 5 = 1/2. 
When po = the change of variables 



a = (uv) 



1/4 



and 



and the new time coordinate dr = adr], allow the field equations (10), (11) and (12) to be recast 
in the simple form 



u = 
2. 



v 



3 
2 



uv = - (V u - k) . 
o 



Solving these equations yields 
u = ci 3 (r - c i4 ) 



and 



u = ^(t - CujVfl - ^ — (r - c 15 ) 
3 3 ci 3 



which on transformation back to the original variables gives 



.2 ^ , J , ^0 2 

as = h < Ci 6 - KT + — T 

r 2 



1 — nr 2 



where c±q = /tci3(c 15 — c^) and all other constants have been absorbed into coordinate redefi- 
nitions. In this case we can transform time coordinates back to proper time t, giving 



L V i — Kr / 



(20) 



As well as this solution there exists a second with 

K 



u = 



V 



and 



V 2 

v = cis + ci 9 r + y r . 



Transforming this solution back gives Milne space. 



4 Perfect fluid solutions 



Perfect fluid dominated solutions to equations (10), (11) and (12) will now be presented. These 
solutions are for spatially flat (k = 0) cosmologies. 
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4.1 d = 1/(37 - 1) 



For the case 5 = 1/(37 — 1) we can integrate the field equations exactly, to find the general 
solution. Notable exceptions are given by 7 = and 7 = 1/3. For 7 = 0, often associated with 
vacuum energy, the corresponding value of 5 is —1, which gives a gravitational theory derived 
from a Lagrangian of the form R° =constant. Clearly this is of little interest. For 7 — » 1/3 the 
value of 5 — > 00, and we have a theory that is not clearly defined. Having excluded these cases 
we continue by changing the time coordinate to 

dr = dr]a z e~ H , 

and considering the variables u and v defined by the transformations 

In v 



(1 = U 2 ( 3 +v / 3 A ) U 2(3— /3A) 



and 



ln-u 



(V3-A) (V3 + A)' 

These definitions force us to exclude A = ±\^3 (corresponding to 7 = —1/3 and 7 = 5/3). 
When A = 2/^/3 — \f?> r ) (corresponding to 5 = 1/(37 ~~ 1)) we can then write the field equations 
(10), (11) and (12) as 



1 1+67 
u = -^7(5 - 37)p ^ 1+37 

v = 

w = ^(5-37)(l + 3 7 ) (V + 



(21) 



37 \ 

p v 1+37 J . 



Solutions to these equations are 



u = 1(5 - 3 7 )(1 + 3 7 ) 2 po(r - c 20 ) 2 (c 21 (r - c 20 ))~^ + (5 3 ^ 1 + 3 ^ V ° ( r - C22 ) 



6 

v = c 2 i(r - c 20 ). 



6c 
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Transforming back to the original variables, and absorbing constants into coordinate redefini- 
tions, gives the metric 



6(1-7) 
5 — 37 



ds 2 = -t TTW ( (1 + 3 7 ) PoT (i+37) + v (t + c 23 )J " 1 dr 2 

2 

+ ((1 + 37)p T< T ^) + V (r + c 23 ) j 5-37 rfx 2 



(22) 



where c 23 = c 2 i(c 20 — c 22 ). This solution is valid for all 7 7^ —1/3, 0, 1/3 or 5/3. 

A second solution to the equations (21) also exists, when Vq < 0. This solution is given by 

1+67 



1 



u = c 24 + c 25 r - -7(5 - 37)p - 

! V Po 



v = 



1+37 

Po 
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Which corresponds to the metric 



6(1-7) 

1 \ 5-37 

ds 2 = - ( (c 24 + c 25 r) - -7(5 - S-y)p T 2 J dr 2 

I \ 5^37" 

+ ( (c 2 4 + c 25 r) - -7(5 - ^)p T 2 I dx 2 



(23) 



Again, this solution is valid for all 7 7^ —1/3, 0, 1/3 or 5/3. 



4.2 * = -(4-37)/(2(7-37)) 



For the theories 5 = — (4 — 37)/(2(7 — 37)) we can also integrate the field equations directly, 
to obtain general solutions for spatially flat cosmologies. We should note that the particular 
cases 7 = 4/3, 7 = 7/3 and 7 = 10/3 are not usefully solved for here. Unfortunately, these 
exceptional cases include the physically interesting case of 7 = 4/3, a fluid of black-body 
radiation. Here 7 = 4/3 corresponds to 5 = 0, which is GR. The case 7-^7/3 corresponds 
to the limit 5 — > 00, and 7 = 10/3 corresponds to 5 = —1. Neither of these values of 5 are of 
physical interest. 

Making the transformation to the variables u and v via the definitions 



and 







1 2 ( -^=v + u 
\/3 



and redefining time as dr] = a 3 dr allows the field equations (10), (11) and (12) to be written as 



u 



(4 - 37) 4(5-37) , 



12 

1 24 (5~37) 

2 



• 2 -2 
U — V 



(4-37) 2 4(5-37) „ (4 -37) 2 Jr - 24 ( 5 -^v 



24(5-37)' 



where we have taken 



A 



24(5 - 37) 
3(2-7) 



V3 (4-3 7 )' 

which corresponds to 5 = —(4 — 37) /(2 (7 — 37)). The case 7 = 5/3 does not give a sensible 
limit, and so we will exclude it from consideration. When 7 7^ 5/3 solutions to these equations 
are given by 



u 



V = 



(4 - 3 7 ) 
4(5 - 37) 
(4 - 3 7 ) 2 
24(5 - 37) 



In 
In 



(5-37)p . 



fir 2 
uc 26 



sin 2 (c 26 (r - c 27 )) 



(5-3 7 )Vb . 



fir 2 
uc 26 



sin 



6c 26 
(4 - 3 7 ) 



[t - c 28 ) 
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The two sin functions here can be transformed to cos functions by suitably redefining the 
constants C27 and C28- These sin or cos functions can then be transformed to sinh or cosh 
function by the transformation C26 — > ic-2%- It is, of course, important that solutions remain 
real for at least some range of r. This restricts which functions should be taken as physically 
interesting (assuming po > 0), so that the form of the corresponding metric depends upon the 
sign of 5 — 37. We will treat the different cases separately below. 



4.2.1 5-37>0 

The condition 5 — 37 > includes the important cases of pressureless dust and vacuum energy 
density, and leads to the metric 



ds 2 = -b 2 (r)dr 2 + a 2 (r)dx 2 



(24) 



where i — 1, 2 or 3. We then have 



1 (4-37) 
ai = sin (5-37) I r — c 29 | s in 6 ( 5 - 3 t) 



(7-3 7 ) (l--y)(4-3 7 ) 

bi = sin (5-37) I r — c 29 | s in 2(5-37) 



6 



(4 - 3 7 ) 
6 



(4 - 3 7 ) 

where constants have been absorbed into coordinate redefinitions, and C29 = C2&{c2j — C2s). 
These sin functions can be transformed to cos functions through redefinitions of the origin of 
the time coordinate r, and the constant C29. A second solution also exists with 

6 



(4-37) 



02 = sinh t 5 - 3 ^) {r — C29} sinh 6 ( 5 - 3 t) 



(7-37) 



(l-7)(4-3 7 ) 



b 2 = sinh < 5 - 3 ^) {r - c 29 } sinh 2(5-37) 



and a third with 



(4-37) 



as = sinh < 5 - 3 t) {r — c 30 } cosh 6(5-37) 

(7-37) _ 

63 = sinh (5-37) | r — c 30 } cosh 2 ( 5 - 3 t) 

These solutions are valid for all 7 < 5/3 and 7^ 4/3. 



(4 - 3 7 ) 
6 

(4 - 3 7 ) 
6 



(4-3 7 ) 

;J -"-I ■'-) f 6 



T 



(4 - 3 7 ) 



4.2.2 5 - 37 < 

The condition 5 — 37 < contains the important case of a scalar field, 7 = 2, and gives the 
metric 



ds 2 = -b 2 Ar)dr 2 + a 2 Ar)dy* 



(25) 



12 



where j = 4 or 5. Here 



1 (4-3 7 ) f 6 

04 = cosh (s- 3 t) {r — C31} sinh 6(5-37) I — -7- 



(4 - 3 7 ) 

(7-37) r , (l-7)(4-3 7 ) f 6 I 

64 = cosh (s- 3 -') {r - c 3 i| sinh 2 ( 5 " 3 t> < — — — r ^ 

1(4-37) J 



and a second solution is 



1 _ (4-37) 

a 5 = cosh < 5 - 3 t) {r — c 32 } cosh 6 ( 5 - 3 t) 

(7-3 7 ) (l-7)(4-3 7 ) C 6 

o 5 = cosh ( 5 ~ 3 t) {r — c 32 } cosh 2 ( 5 - 3 x> — -r 



(4 - 3 7 ) 
6 

(4 - 3 7 ) 



These solutions are valid for all 7 > 5/3 and 7^ 7/3 or 10/3. 



5 Analysis of vacuum cosmologies 

In this section we will perform an analysis of the solutions found in section 3. The special cases 
of 5 = —1/2 and 5=1/4 will be investigated in the appendix. These solutions contain a number 
of constants, as a result of integrating the field equations. The physical significance of these 
constants varies: In some cases they can be absorbed into a rescaling of coordinates, in others 
they cannot and must be specified by initial conditions 2 . If the latter is the case, then these 
constants are physically meaningful quantities and their value is important for the evolution of 
the space-time. We will describe the effect of taking different values for these constants in the 
analysis below. 

5.1 Spatially flat cosmologies 

We begin with the spatially flat solution (19). In the case en = the solution (19) can be 
written as 

ds = -e (is-m2s+i) T dT 2 + eW-v dx. 

2,5(1+2,5) 

= -dt 2 + t V=V dx 2 , 

which is simply the power-law solution (8). It can be seen directly from (19) that for all en 
this simple solution is the attractor as r — > 00. 

When en = the power-law solution above describes the evolution of the universe all 
the way back to the initial singularity, but for en 7^ more general behaviour occurs as the 
singularity is approached. The form of this early-time behaviour depends upon the sign of c n 
(recall that the magnitude of cn can be absorbed into coordinate redefinitions). 



2 This is not a new phenomenon, and can be shown to be the case in, for example, the general Friedmann 
solutions of Brans-Dicke theory [30] . 
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For cu > a power series expansion of the metric gives 



ds 2 ~ -T^-^dr 2 + TT&^dx 2 
= -dt 2 + trfx 2 

when t -\ log(cn) <C 1. This limit is the same as the power- law exact solution (9). 
For en < a similar expansion about |log(— c n ) gives 




, 9 / 5 (i+2g) \ 2 „ 



This solution is clearly quite different from the power-law limits that have been found so far. For 
<5>l/4or<0it corresponds to a non-zero minimum of expansion, or bounce. For < 5 < 1/4 
it is a maximum of expansion. It should be noted that when — 1/2<5<0 the power of t in the 
scale factor diverges as t — > 0. In this case there is no staticity, and the scale-factor diverges as 
t — > 0. Bounces occur here without any violation of energy conditions, as there are no matter 
fields present to commit such violations. The bounce is entirely due to the vacuum dynamics 
of the theory. 

We have now shown both the early and late-time behaviour of the solution (19). At late- 
times all solutions approach the power- law exact particular solution (8). At early-times the 
evolution of the universe is qualitatively different depending upon the value of the constant 
c\\. For cu = the solution (8) is valid all the way back to the initial singularity. For 
Cu > the general solution has an early period of expansion of the form of (9), in the vicinity 
of the singularity. The expansion then evolves into the late-time attractor. When c n < 
there is a non-zero minimum of expansion (or maximum if < 5 < 1/4, or divergence if 
— 1/2 < 5 < 0). The universe then evolves from this static past towards the late-time attractor. 
Some representative evolutions of the scale-factor, in terms of the proper time coordinate t, are 
shown in figure 2. 

We will now investigate the solution (18). Firstly, expanding around the point r = we get 
the approximation 



when r < 1. This is remarkably similar to the form of the previous metric, (19), in the vicinity 
of its minimum of expansion. A noticeable difference is the sign before the second term in the 
brackets. This change in sign shows that when the previous solution was a minimum (5 > 1/4 
or < —1/2) this solution is a maximum of expansion. Correspondingly, when the previous 
solution was a maximum (0 < 5 < 1/4), this solution is a minimum. Again, there is divergence 
when -1/2 < 5 < 0. 
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Figure 2: The evolution of the scale factor of solution (19) in terms of proper time, t. The value 
of S chosen here is 1/2. The solid line corresponds to cn = 0, or the power-law exact solution 
(8). The dotted line corresponds to c n > and shows an early period of t 1 ! 2 expansion, that 
is attracted towards (8) at later times. The dashed line corresponds to c n < and shows a 
non-zero minimum of expansion, and subsequent evolution towards (8). 

We will now investigate the form of (18) around r = ±7r/2. Performing a power series 
expansion about either of these points gives 

ds ~ -rW^dr 2 - r^^dx. 
= -dt 2 + tdx 2 , 

when r =F 7r/2 <C 1. In the vicinity of both of these points the evolution of (19) is therefore of 
the form of the power-law exact solution (9). 

We have shown that the evolution of the solution (18) depends upon the sign of 5. When 
5>l/4or<— 1/2 the solution, in terms of proper time t, has a maximum of expansion with its 
early and late-time evolution going as frs. When < 5 < 1/4 there is a minimum of expansion, 
and late and early-time expansion again goes as t l l 2 . For — 1/2<5<0 the solution diverges 
in the region which is static for all other values of 5. The form of some representative solutions 
are shown in figure 3. The behaviour of these solutions are markedly different from those found 
by Russo [28]. Although a late-time power-law attractor can still be seen to exist in one of the 
solutions, very little else is comparable. 

5.2 Spatially curved cosmologies 

For spatially curved vacuum universes we found the general solutions (20) and Milne space, 
when 5 = 1/2. The existence of Milne space as a solution should not be surprising, as any Ricci 
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Figure 3: The evolution of the scale factor of solution (18) in terms of proper time, t. The solid 
line corresponds to 5 = 1, the dashed line to 5 = 1/8 and the dotted line to 5 = —2. 

flat solution of GR in a vacuum is also a solution of the field equations (2) in vacuum (at least 
when 5 > 0). We will now investigate the form of the solution (20). 

When k = this solution reduces to 

ds 2 = -dt 2 + t 4 dx 2 , 

which corresponds to the power-law particular solution (8), when 5 = 1/2. When K^Owe see 
that this solution exhibits new behaviour, dependent on the value of n and which branch of the 
± sign is chosen. 

When k > and cn > 0, then in the vicinity of t = this solution exhibits a maximum 
of expansion. Subsequent evolution of the scale-factor, as t — > ±oo, depends upon the chosen 
branch. Taking the positive branch we see that the universe can start to expand again, as 
\t\ increases. The existence of this expansion phase depends upon the magnitude of c 17 . If 
Cn ^ k 2 /4 then the effect of the second term in the scale-factor is too great and the universe 
collapses to a singularity, before the t 4 term can become dominant. When cyj > k, 2 /4 the 
t A term dominates the late-time evolution and the power-law solution (8) is approached, as 
\t\ — > oo. When the negative branch is taken then the maximum of expansion at t — always 
leads to collapse to a singularity, as t either increases or decreases. 

When k < and c 17 = then there exists an initial singularity, with the evolution of the 
scale factor in its immediate vicinity expanding proportionally to t. The late-time evolution of 
such solutions will be attracted to either the power-law solution (8), if the positive branch is 
chosen, or to eventual collapse to singularity, if the negative branch is chosen. The case cn > 
does not have a singularity at t — 0, but instead has a minimum of expansion, or bounce. The 
evolution away from this minimum proceeds as in the cyj = case. The cn < case again 
features an initial singularity and displays approximately the same evolution as when cn = 
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(a notable deviation is that the period of a oc t expansion near the singularity is absent if cyj 
is sufficiently negative). 

We have shown that when k 7^ there exists a minimum of expansion when Cyj > 0, and 
an initial singularity when cyj ^ 0. Initial evolution away from an early minimum proceeds 
as a oc t, either expanding or collapsing, depending on the sign of k. If there is an initial 
singularity, then evolution away from it proceeds as a oc t, unless cyj is sufficiently negative (in 
which case the t 4 term dominates right from the beginning). The late-time evolution of these 
solutions is dominated by the t 4 term in (20). The negative branch of this term causes collapse 
to a singularity, and the positive branch leads to asymptotic expansion and approach towards 
the power-law solution (8) (unless the term proportional to k causes collapse to singularity 
before the t 4 becomes dominant). The form of some representative solutions are shown in 
figure 4. 



a 




Figure 4: The evolution of the scale factor of solution (20) in terms of proper time, t. The 
solid lines correspond to universes with k > 0, dashed lines to those with k < and the dotted 
line to k — 0. Those solutions starting from the origin have Cyj = 0, and those starting from 
a non-zero minimum have Cyj > 0. All solutions collapsing to singularity correspond to the 
negative branch of t 4 , and all solutions expanding at late-times to the positive branch. See the 
main body of text for an explanation of these different behaviours. 



6 Analysis of perfect fluid cosmologies 

We will now analyse the evolution of universes described by the perfect fluid solutions (22), 
(23), (24) and (25). Again, these solutions contain a number of constants, some of which 
can be absorbed into coordinate redefinitions and some of which cannot. We will explain the 
significance of such constants below. 
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6.1 d = 1/(37 - 1) 



When 5 = 1/(37 — 1) we found from direct integration of the field equations the general solutions 
(22) and (23). We will first consider (22). 

The form of the scale factor in this solution depends upon which of the two terms in the 
brackets is dominant. When the first term dominates the metric is approximately given by 

„ 6(l+4 7 -3 7 2 ) 2 

ds ~ — t (5-37x1+37) dr + t (5-37x1+37) dx. 
= -dt 2 + tdx 2 , 

which corresponds to the power-law particular solution (9). When the second term in the 
brackets dominates we have 

6(l+7 7 -6 7 2 ) 2 „ 

ds ~ — t (5-37x1+37) + r (5-37x1+37) rfx 

2(l+3 7 ) „ 
= -dt 2 + t (2-37X1-37) rfx 2 , 

which corresponds to the vacuum dominated, power-law particular solution (8). It now remains 
to investigate the conditions under which each of these approximations dominates, at early and 
late times. 

When C23 = 0, and 7 > 0, then the first term always dominates as r — > 0, and the universe 
evolves like a ~ t?. As r — > 00 the second term dominates, and the evolution of the universe is 
attracted towards the form of (8). 

When c 2 3 7^ 0, and 7 > 0, then this has the effect of offsetting the origin of the second term 
relative to the first. The second term will then be the dominant one for most of the range of r, 
with a contribution expected from the first term in the vicinity of r = — C23. 

Phantom fluids, with 7 < 0, have the same behaviour as outlined above, but with the 
dominance of the two terms reversed in each case. 

We will now consider the solution (23). The evolution of this space-time again depends on 
which of the two terms in the brackets is dominant. If the first term dominates then the metric 
look like 

, o 6(1-7) „ 2 

ds 2 ~ -t W^Jdr 2 + r^^Tldx 
= -dt 2 + tdx 2 , 

which is again the power-law evolution (9). We will call this vacuum domination, as it cor- 
responds to evolution being dominated by a term independent of p . If the second term in 
brackets dominates then we now have 

12(1 — 1) 4 

ds 2 ~ -t~ (5- 3 7) dr 2 + rT^ldx 2 
= -dt 2 + tV^h)dx 2 , 

which corresponds to the matter dominated power-law expansion (7). We must now investigate 
the conditions under which each of these approximations dominates. This will depend upon 
the value of the constants c 2 4 and c 25 . 
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When C24 = C25 = we see that the matter dominated power-law expansion (7) holds all 
the way back to the initial singularity. If C24 = and C25 7^ 0, then there is an initial period of 
vacuum dominated expansion of the form a oc £2 and a late-time evolution towards the matter 
dominated power-law solution (7). If C24 7^ then this has the effect of offsetting the origin 
of first term in the brackets, relative to the second. This can result in a non-zero minimum 
of expansion (depending on the signs of C24 and C25), as in this case the first term can both 
dominate and be non-zero at r = 0. As \r\ — > 00 the matter dominated power-law solution (7) 
is then approached. 



6 .2 5 = -(4-3 7 )/(2(7-3 7 )) 

We will now investigate the solutions (24) and (25), obtained by integrating the field equations 
in the presence of a perfect fluid when 5 = —(4 — 37) /(2 (7 — 37)). 

6.2.1 5-37>0 

We will first consider the solution {ai,&i}. Expanding this solution around the points r = nn 
and r = c 2 g + nn, where n is an integer, gives the evolution of the scale factor in their vicinity. 
To first order, an expansion about t = nn results in 

(4-3 7 )(l-7) (4-3 7 ) „ 

ds ~ -r ( 5 - 3 7) dr l + r (b-3-t) dx z 

2(4-3 7 ) , 
= ~dt 2 + £3(7-37X7-2) rfx 2 , 

which corresponds to the vacuum dominated power-law expansion, (8). A similar expansion 
about r = c 2 g + nn gives 

2(7 — 37) 2 

ds 2 ~ -(r - C2 9 y~^^dT 2 + (r - c 29 )"(^^c/x 2 
= -dt 2 + tdx 2 , 

which is the power-law evolution described by (9). Whether or not these points correspond to 
a curvature singularity depends upon the value of 7, and can be read off from the Ricci scalar 

6(10 - 37) 2(7-37) f , (2- 7 )(7-37) f 6 

R = i — ttt 21 — r sin v>--A~o { r - c 29 } sin ( 5 - 3 ^) ' 



(5-3 7 )(4-3 7 ) " ( ' zaJ ^ l(4-3 7 ) 

We see that when 7 ^ 7/3 the Ricci scalar remains finite for all r, whilst for 2 < 7 < 7/3 or 
7 < 5/3 the Ricci scalar diverges to infinity every time r = (4 — 37)n7r/6. The evolution of 
a(t) in the vicinity of these singularities goes like (8). Similarly, when 5/3 < 7 < 7/3 there are 
curvature singularities at r = c 2 g + nn. In the vicinity of these singularities a(t) ~ £2. These 
latter singularities correspond to a(t) — > 0, whereas the former correspond to the divergence 
a(t) — > 00 (except in the range 4/3 < 7 < 5/3, in which case a(t) — > 0). 

An instructive special case to consider is that of pressureless dust, 7 = 0. In this case the 
coordinate transformation £ = cot(r — c 2 g) allows the metric to be recast as 



ds 2 = —dt 2 -\ j — c£x 2 

sin^tcag + cot- 1 ^})} 
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It can now be seen directly that a(t) — > t 2 as t — > oo (or r — > c 29 + rar), and diverges to oo as 
(8) when t — > — cot{c2g — mr/6} (or r — > mr/6). 

We will now consider the solution {02, 62}- In the limit r — > ±00 this solution approaches 

ds 2 ~ -e^Is^Trfr 2 + e'^^Idx 2 
= -dt 2 + tdx 2 , 

which is the power-law solution (9). When r^Owe see that we have the same limit as in the 
previous solution, where the power-law evolution (8) is approached. Similarly, as r — > c 2 g we 
approach the power-law evolution a ~ t?, (9). The conditions for whether or not these limits 
correspond to singularities, and if these singularities correspond to a(t) — > or a(t) — > 00, have 
exactly the same dependence on 7 as with the previous solution. 

The pressureless dust case is again instructive, and can be written in terms of proper time, 
by making the transformation t = coth{c29 — r}, as 



ds 2 = —dt 2 H j dx 2 . 

sinh6{6(c 29 - coth -1 -^})} 

This solution displays evolution as t — > 00, and divergence to 00 as t — > coth{c2g}. 

The solution {03, 63} will now be considered. The late-time evolution of this solution is the 
same as in the previous case. The evolution about r = c is also the same. Now, however, we no 
longer have the singular behaviour about r = that previously existed. Instead the scale-factor 
evolves as a ~ a + t2 in the vicinity of this point, where a is some constant. The dust solution 
in this case is then 

Jfi - 1 

ds 2 = -dt 2 + j dx 2 , 

cosh6{6(c 30 - coth 1 {t})} 

where t = coth{c3o — r}. This metric has a scale factor that no longer diverges at any finite t. 



6.2.2 5 - 37 < 

The late-time evolution of these two solutions is of the form a ~ t^, as with the previous two 
solutions. Now, the form of solution {04, 64} about r = displays the same divergences as {02, 
b 2 }. Evolution towards the point r = c is non-singular, and goes as a ~ a + t*. The special 
case of pressureless dust is given in this case, in terms of proper time t, by 



ds 2 = —dt 2 H j dx 2 , 

sinh6{6(c 3 i + tanh~ 1 {t})} 

where t = tanh{r — c 31 }. The scale-factor here can be seen to vanish at t — ±1 (corresponding 
to t — > ±00), and to diverge to 00 at t = tanh{— c 3 i}. The evolution of the scale factor follows 
the late-time attractor towards the points at t — ±1. 

We will now consider the remaining solution, {a 5 , 65}. This solution behaves as the solution 
{a 3 , 63} in the vicinity of r = 0, and as {a 4 , 64} in the vicinity of r = c. That is, about both of 
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these points the evolution of the scale factor goes as a ~ + and is non-singular. In fact, 
at no finite r does this solution become singular. The pressureless dust solution is given here, 
in terms of proper time, as 



ds 2 = -dt 2 H ; rfx 2 , 

cosh6{6(c 32 + tanh -1 {£})} 

where t = tanh{r — c 32 }. This solution describes a space-time that evolves from a = at 
t — — 1 (corresponding to r = —00) to a — at t — 1 (corresponding to r = 00). The evolution 
between these two end points is at all times finite. 



7 Discussion 



We have investigated the homogeneous and isotropic cosmological solutions of R n theories. 
It has been shown that in a number of cases the field equations can be integrated directly, 
allowing the general behaviour of these models to be found. For spatially flat vacuum universes 
the solutions for any n can be found. For spatially curved vacuum universes, and flat perfect 
fluid universes, the general solutions can be obtained for various particular values of n. These 
solutions were given explicitly in sections 3 and 4. It is interesting to note that unlike the 
general relativistic cosmologies, vacuum solutions can be either open, closed or flat. It was 
also found that there is generically not one unique solution for any given n, matter content 
and topology, as there usually is in GR. This was found previously in [17, 21], where a phase 
plane analysis showed an invariant sub-manifold in the phase space of solutions which could 
not be crossed by any trajectory. This result is reiterated here, where we find multiple ways 
to integrate the field equations, yielding multiple solutions. It was also shown in [17, 21] that 
the form of the general solution changes with n: Critical points in the phase space change their 
attractor nature as n is varied. Again, we confirm this result by finding explicit solutions for 
cosmologies with the same topology and matter content, but different n. These solutions have 
clearly different forms. In sections 5 and 6 we performed a brief analysis of the vacuum and 
perfect fluid solutions, respectively. The evolution of these universes, in terms of proper time, 
were found in different regimes. The nature of the initial singularity was investigated, or its 
lack thereof, as well as evolution at late-times. In summary, it was found that great variety is 
present in the cosmologies of these theories. Behaviour at both late and early times can vary 
wildly, depending on n, and even between different solutions with the same n, matter content 
and topology. Solutions can manifest initially singular or non-singular behaviour, and late-time 
evolution can either lead to a crunch or to monotonic expansion. 
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A Vacuum S = —1/2 cosmologies 

We will now describe the evolution of universes described by the solution (16). Taking the 
positive branch of (16) we see that as r — > oo we have dr <C dt, so that in terms of proper 
time the universe evolves very slowly at late-times. However, for the negative branch of (16) 
we see that dr ^> dt in this limit. For the negative branch, therefore, we should expect 
increasingly rapid evolution at late times, and for the positive branch an evolution towards an 
asymptotic steady state. A power series expansion of (16) shows that the early-time behaviour 
of both branches goes like the power-law particular solution (9). Figure 5 shows graphically 
the evolution of the two branches of this solution. 



a 




Figure 5: The evolution of the scale factor of solution (16) in terms of proper time, t. The solid 
lines correspond to the positive branch, and the dashed line to the negative branch. 



B Vacuum S = 1/4 cosmologies 

We will now describe the evolution of universes described by the solution (17). A power series 
expansion of this solution shows that the scale factors of both branches approach a stationary 
state, in terms proper time t, as r — > 0. We can then see that as r increases the positive 
branch expands out of this stationary state, and the negative branch collapses out of it. The 
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positive branch therefore corresponds to an expanding universe evolving out of a bounce, and 
the negative universe to an initially expanding universe with a maximum of expansion and 
subsequent collapse. Figure 6 shows the evolution of the two branches of this solution. 




t 



Figure 6: The evolution of the scale factor of solution (17) in terms of proper time, t. The solid 
lines correspond to the positive branch, and the dashed line to the negative branch. 
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